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Disoriented Chiral Condensates (DCCs), the correlated space-time regions where the chi- 
ral order parameter of QCD is chirally rotated from its orientation of the true vacuum state 
in the isospin space, might be formed in high energy hadronic or nuclear collisions. The 
search of such a DCC would provide a probe to the understanding of the chiral structure of 
the QCD vacuum and/or the chiral phase transition of strong interactions at high temper- 
atures. Various experimental searches of this phenomenon include cosmic ray experiments, 
the MiniMAX experiment in proton-antiproton collisions at the Fermilab Tevatron, or event- 
by-event analysis of Pb-Pb collisions by the WA98 and NA49 collaborations at the CERN 
SPS. Although no clear experimental evidence for DCC formation has been reported so far, 
the search of this phenomenon is still part of the missions in the heavy-ions physics program 
in present or planned experiments at RHIC and LHC [2| . 

In relativistic heavy ion collisions, the large energy deposit in the central rapidity region 
leaves behind a hot and dense plasma at a temperature above 200 MeV where the chiral 
symmetry is restored. Then, the plasma cools down in a rapid hydrodynamic expansion 
accompanied by the chiral phase transition. The subsequent out of equilibrium evolution 
triggers an exponential growth of long-wavelength fluctuations via the spinoidal instabil- 
ities yHsL resulting in the formation of disoriented chiral condensates (DCCs) (see also 
Refs. p, lZI). The decay of such DCCs to the true QCD vacuum is expected to radiate 
copiously soft pions that could be a potential observational signature of the chiral phase 
transition. However, the emitted pions will strongly interact with the background of other 
hadronic matter so that the hadronic signal may be severely masked. In contrast, electro- 
magnetic probe such as photon and lepton with longer mean free path in the medium serves 
as a good candidate and can reveal more detailed non-equilibrium information on the DCCs 
with minimal distortion js- lol . 

It has been proposed that anomalous radiation of low-momentum photon and/or low- 



mass dileptons can be produced from a DCC. In particular, in Ref. jlOj, Boyanovsky et al. 
have extensively studied the photon production from the DCC with an nonvanishing expec- 
tation value of the neutral pion through the Ua(1) anomalous vertex. Later, the authors 
in Ref. [lH have taken into account another dominant contribution that also involves the 
dynamics of 7r° due to the decay of the vector meson through the electromagnetic vertex. 
They have found that for large initial amplitudes of the vr" mean field the photon production 
is enhanced by parametric amplification, resulting in a distinct energy distribution of the 
produced photons. However, the above two works have ignored the hydrodynamical expan- 
sion and adopted the simple "quench" approximation for the chiral phase transition. This 
quenched phase transition has been widely used in the study of non-equilibrium phenomena 
of DCCs S&Sl. Based upon the Bjorken's scenario, in ultra-relativistic heavy-ion colli- 
sions, an approximate Lorentz boost invariant particle flow along the longitudinal direction 
might be created in the central region [l2[. At late times following the heavy nuclei colli- 
sions, a transverse flow can be generated due to the multi-scattering between the produced 
particles, as such the expansion becomes three dimensional 13||. In Ref. the authors 



have considered the effect of the hydrodynamical expansion of the plasma to the production 
of photons from the non-equilibrium relaxation of a DCC. It is found that the expansion 
smoothes out the resonances in the process of parametric amplification such that the non- 
equilibrium photons are dominant to the thermal photons over the range 0.2 — 2 GeV. This 
work is to assume a spherical boost invariant hydrodynamical flow which for small values of 
the rapidity is conformally fiat, thus greatly simplifying the treatment of the photon field in 
the expanding spacetime. In this work, we attempt to consider the nonequilibrium photon 
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production from a DCC in a more realistic longitudinal boost invariant hydrodynamical 
expansion, which certainly requires a quantization scheme for electromagnetic fields under 
this anisotropic expanding background. 

The study of the quantized electromagnetic fields interacting with the gravitation fields or 
propagating in the curved spacetime is of interest in many different aspects, one of which is 
to consider t he p hoton production in the early universe with its possible effect on the initial 
anisotropics [15[. Thus a consistent quantization scheme of electromagnetic theory under an 
anisotropic expanding background is required. In Ref. [l6[, quantization of electromagnetic 
fields in a diagonal Bianchi type I metric is proposed using the Fourier mode expansion of 
the fields developed in Ref. 131 . Applying the WKB particle concept in asymptotic in and 



out regions of spacetime where the background expansion is adiabatically slow gives the nor- 
malization of the mode functions. Then particle production of the free photons is considered 
by performing the standard Bogoliubov transformation between annihilation and creation 
operators in in and out regions. Here we instead try to provide an alternative quantization 
scheme directly on electromagnetic potentials. An application of this quantization method 
to interactive photon fields is straightforward, in particular when some other charged matter 
fields are also present with the interaction term depending upon electromagnetic potentials. 

After introducing our model and summarizing the dynamical equations for describing 
the evolution of the expectation value of the fields within the DCCs, in this Letter we will 
propose the quantization method on quantizing electromagnetic potentials in an anisotropic 
expanding background and discuss its application to the interactive photons. The formalism 
for considering photon production will be developed. Summary and discussions on our future 
investigation will be presented at the end. 

The longitudinal expansion, say along the z axis, of the boost invariant hydrodynamical 
flow can be described by the proper time r and the space-time rapidity rj, defined as 

.4.n(^), (1) 

where (t, x) = (t,x_|_,2;) = (t,x,y,z) are the coordinates in the laboratory, 

t = rcosh?7, 2; = rsinh?7. (2) 

The ranges of these coordinates are set to be < r < oo and < ?7 < oo, restricted to the 
forward light cone. The Minkowski line element is then given by 

ds'^ = dt^ - dy? 

= dr^ - rfxl - r^dr/^ (3) 

Eq. ([3]) is the Kasner metric that serves as the background comoving frame under which 
non-equilibrium photons are to be produced from a DCC domain. 

The dynamics of the DCC as well as the photon production from a DCC in a generally 
expanding background can be described by the phenomenological action given by 

S = J d^x^ {L, + La + L^oa) , (4) 



where 



L. = -^g^'^d^^ ■ 9,$ + ■ $ - A ($ • $) +ha, (5) 
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La 

LnOA 

rO 



IT 



mfmr 



(6) 

(7) 



$ = (cT, vr'^, vf) is an 0{N + 1) vector of scalar fields with tt = (vr^, tt^, 7r^~^) representing 
the N — 1 pions. The large- approximation will be implemented later to account for the 
strong coupling effects among the a and all other vr fields. The phenomenological parameters 
in the effective Lagrangian above can be determined by the low-energy pion physics. F, 



dfj_Au — dyA^ is the electromagnetic field tensor. The dominant couplings between the photon 
and the neutron pion are given by the anomalous Ua{^) interaction as well as the decay of 
the vector meson V^. is identified as the u meson with its coupling Xy that is obtained 



from the u — > 7r°7 decay width [18|. The signature of the metric we choose is ( — h ++) 
where g = — detf^f^j^] = r^. In particular, since e°'^'^'^ is a tensor density of weight —1, each 



of this pseudo tensor is multiplied by the factor in a general coordinate system [l9 



The above effective interactions on the photon have been studied by us and are obtained 
from the perturbative theory without involving in-medium modifications. The in-medium 
effects will enter only through the non-equilibrium fluctuations for the pion flelds due to 
their strong couplings . 

Within the DCC, we assume that the a and the vr'' acquire the respective expectation 
values that both are a function of the comoving time r only in such a boost invariant 
background. Here we flrst shift a and by their expectation values with respect to an 
initial non-equilibrium states: 



7r°(x,r) = 
with the tadpole conditions: 

(X(x,r)) 



+x(x,r), 
C(r) + ^(x,r), 



(a(x,r)) = 
(7r°(x,r)) 



C(r) 



0, (^/'(x,r)) = 0, (7r(x,r)) = 0. 



(8) 
(9) 

(10) 



Using the Schwinger-Keldysh closed-time-path formulation of non-equilibrium quantum fleld 
theory allows us to derive the evolution equation of the non-equilibrium expectation values 
as well as the correlation functions of quantum flelds. This technique has successfully been 
employed elsewhere within many different contexts and we refer the readers to the litera- 
ture for details 20|. To self-consistently incorporate quantum fluctuation effects from the 



strong a — IT interactions, we then implement the large- approximation that provides a 
nonperturbative resummation scheme below (Eqs. (2.7)-(2.9) of Ref. 10|). Using the tad- 



pole conditions (fTOj) . the large- equations of motion for the mean flelds that also involve 
the perturbatively electromagnetic corrections can be obtained as follows: 



T 



-dr 



^ + 4A02(r) + 4AC^(r)+4A(7r2)(r) 



2 

^2 



^ + 4A02(r) + 4AC'(r)+4A(7r2)(r) 



0(r) 
C{r) 



/i = 0, 



327r2/. 



(FF)(r) + 



Amlmy 



gas{F^uFap){r) 



r-'Xe'"''"'e"''^'g^s{d,F^uFap){T) = 0. 



Ami 



ml 



(11) 
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where the dot means d/dr. The backreaction from the quantum fluctuations of the pions 
can be expressed by their Fourier mode functions Ukir) defined as: 



vr X, T 



d^k 



dkUJr) e 



ikx 



h.c. 



(12) 



where k-x = k_|_x_|_ + k^?7 and d^k = d^k^dk^^. is a destruction operator for the mode with 
frequency a;^^(ro) to be determined from solving the time independent gap equation given 
by the initial states. The mode equations for the pions can be read off from the quadratic 
part of the effective Lagrangian as 



M: 



Uk{T) = 0, 



(13) 



where k"^ = g^uk^k'^ = k\ + ^^/t^- The expectation values with respective to the initial 
states are given by 



(7f2)(r) = (iV-l) 



d^k 



2(27r)3a;.,(ro) 



coth 



2T,- 



(14) 



where (7r^)(r) will be self-consistently determined by the above equation with an initial tem- 
perature set to Tj presumably above the critical temperature of the chiral phase transition. 
Thus, the corresponding the gap equation for the initial adiabatic modes is 



ki' 



Ml 
2 



+ 4A02(ro) + 4AC'(ro) + 4A(7f2)(ro) 



(15) 



It must be noticed that this is an effective field theory with an ultraviolet momentum 
cutoff of the order of A my. Without considering the perturbative backreaction effects 
from the photons, the above equations have been extensively studied in Ref. aiming at 
finding the phase space of initial conditions that may lead to instabilities on long wavelength 
modes during the subsequent nonequilibrium chiral phase transition. If so, then the DCCs 
may be formed, and thus produce an anomalous transverse distribution of secondary pions 
when compared to a more conventional boost invariant hydrodynamic flow in local thermal 
equilibrium. Here these background solutions are of importance to provide a nonequilibrium 
environment from which the photons are to be produced. 

The relevant Lagrangian densities for describing the dynamics of photon production under 
the time dependent expectation value of the neutral pion are given by La and L^^ in 
particular where the pion field in L^o^ is replaced by its mean value C,. The underlying gauge 
symmetry embedded in electromagnetic theory reflects the fact that not all gauge potential 
are of physical relevance. Thus we can avoid potential gauge ambiguities by eliminating all 
possible redundant degrees of freedom with a choice of the gauge fixing. 

The canonical momenta, 11^ = (11^,11*), conjugate to the potential fields, = {Ar,Ai), 
are defined by 



OA-r- 



5Ai 



(16) 
(17) 
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However, since A^- is not a dynamical variable as its conjugate momentum n"^ vanishes, we 
can set A^- = 0, a. temporal gauge. Then the further constraint equation can be obtained by 
taking the variation of the Lagrangian density with respect to A^ that gives the Gauss law 
diW = 0. Thus, the choice of the gauge fixing can be such as 

A, = 0; g'^d,A,=0, (18) 

the counterpart of the Coulomb gauge in a fiat spacetime. The gauge conditions leave us 
two transverse vector potentials. It proves more convenient to parametrize the polarization 
vectors in a longitudinal expanding background as follows: 

They obey 

^(1) _ ^(1) ^(2) _ _^{2) Ja) I _ J_ Imn Jb) (c) /^„x 

^-ki ~ ^ki ' ^-ki— ^ki ' ^k ~ y/^ ^km^kn' K'^^J 

where (a, 6, c) = (1,2,3) and its cyclic order. Notice that 

so the longitudinal expansion might lead to the change of the direction of the polariza- 
tion vector e^\^ but keeps the other two polarization vectors e^] and remaining their 
directions. Thus, the mode expansion of the vector potentials must involve not only the 
state with polarization vector e\^l, but also the mixed state with the respective polarization 
vectors and ejf], given by 



h, 1 

. , , { hk (V^u.(r)^i + V^3;^'(r)41) + &2kV^2.(r)eS] + h.c.} . 



d^k 

(22) 

The creation and annihilation operators satisfy the canonical commutation relations. The 
mode function Vsi^^r) certainly is not an independent function and its evolution depends on 
the choice of the gauge that gives 

^3.W + ^^3.(r) = ^Mr). (23) 

It is then expected that the physical degrees of freedom are those two transverse modes with 
the polarization vectors obeying the transversality condition. 



a=l,2 



E4141=^.-^ = ^.(k). (24) 



It will be seen later that the equation fl23l) obtained from the gauge fixing then becomes 
crucial in the sense that the expression of the gauge invariant quantity of interest, for example 
electric and magnetic fields, have no V^k dependence. 
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Then the mode expansion of the conjugate momenta can be obtained as 



(25) 



where the electric and magnetic fields can be expressed in terms of two transverse degrees 
of freedom as follows: 



-(2)i 
-k 



e^'^'' + h.c 



S^(x,r) 



2 

dPk 



./g{ik 



{2)i 
k 



k 



^jkx 



+ h.C.} 



SO that the fields always point to the direction perpendicular to the direction of their prop- 
agation k. 

The commutation relation between the vector potentials and their conjugate momenta 
can be straightforwardly computed as follows: 



+V9 



-ihPUk) 



2k{ro) 



Vikir) 



kH 



k' 



kH 



lk\ 



^k 

-(1) i 



(1) i 



+ 



CikV^k{r)e'^ 



{2)i 



+ 



(2)i 



(3) 



,(2J 



where the Wranskian conditions of the mode functions Via:(t) and V2A;(r) and their time 
derivatives Vifc(r) and V2k{T) are given by 



2k{To) 

1 

2k{To) 



Vik{r)VUr) - VUrmkir] 
V2kir)V,liT) - V*,iT)V,,iT] 



—ih . 



(29) 
(30) 



The Wranskian equations are to fix the normalization conditions of the mode functions. The 
first term of the commutation relation above (128|1 is expected whereas the additional terms 
with the V3fc(r) dependence are present to satisfy the gauge choice, that is, [W^ , y A_y_j ] = 
obtained by means of its equation (123|) . Notice that since we are interested in the dynamics 
of photon production from a DCC that does not necessarily involve the intermediate photon 
states, it will be seen later that the expression of the number density of the produced photons 
depends only on the dynamics of the mode functions respectively, Vife(r) and V2fc(r). 

Although the quantization scheme can be formulated in the way that all redundant de- 
grees of freedom can be eliminated consistently by choosing a proper gauge fixing, there 
might exist a subtlety on finding a suitable mode expansion to further diagonalize La- 
grangian density of the interactive photon fields. This subtlety will become more clear later 
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as the above quantization scheme is apphed to free photons first and to interactive photons 
of interest later. Besides, we will mention the known alternative quantization scheme in 
Ref. 16| that is proposed for investigating the dynamics of photon production under the 
more general time dependent anisotropic spacetime, and comment on the advantage of our 
formulation. 

The Lagrangian density for free electromagnetic fields La can be obtained by imposing 
the gauge conditions (fTg|) as 



1 



1 



^lAT^^'A 



T ) 



(31) 



where the equations of motion for the Fourier modes of the vector potentials A^-k are ob- 
tained as follows: 

dlA^T^ir) + -d^A'^^ir) + k'A'^^ir) = . (32) 
r 

The mode functions for the linear polarization states are thus 



2k] 



r 








(33) 
(34) 



where the Wranskian conditions fl2^ and fl5Ul) are satisfied. In Ref. [16[, the quantization 
scheme is based upon the method developed in Ref. [l7| where the mode functions S'(±)fc to 
be defined later are used in the Fourier expansion of the electric and magnetic fields. The 
mode functions S(^±)k can be related to the mode functions V,*), and as 



= 2 



k^T 



k 



(35) 



Then, using Eqs. (!33|) and (!34|) . S'(±)fc can also be expressed as 



S: 



(±)fc 



kv,l±v,l 



2k 



(36) 



and the evolution equations of S'(-|-)fc are found remarkably simple, that is. 



S{±)k{r) + ^5'(±)fc(r) + A;^S'(±)fe(r) = , 



(37) 



obtained in Ref. [17|. The corresponding Wranskian conditions of the mode functions S(^±)k 
can also be found. Here we propose the quantization scheme on the electromagnetic poten- 
tials directly. One immediate merit is that the normalization condition of the mode functions 
that will be used later to compute the number density of produced photons can be straight- 
forwardly determined. It is particularly useful to set up the initial-valued problem of the 
nonequilibrium system where the initial mode functions can be specified consistently sub- 
ject to the Wranskian conditions. However, this huge simplification works only for the free 
photons. Can neither Vi2k 'S'(±)fc further diagonalize the Lagrangian for the interactive 
photon fields where the appropriate perturbation expansion will be employed later. 
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As for the interactive photon fields under consideration, the corresponding Lagrangian 
density can be written in terms of the transverse degrees of freedom as 



La + hA = Iat^A't - i f 1 - 4^ C'l diArid'j^T ' j^C^e'^'' Ati O^Atu ■ (38) 
2 2 V mlmi, J An^f^ ^ 

lYiom. observing the mode equations fl33|) and f lM|) . we add and subtract the term 

Li-2^Tuy-^,--^,jA^,, (39) 

in the Lagrangian density La + Lc_a- The subtracted term will modify the kinetic energy of 
the transverse modes of the vector potentials with the polarization vector e^^\ so that the 
expansion modes in terms of the circularly polarized states to be defined later can diagonalize 
the unperturbative Lagrangian density given by La + L,^a ~ Lj whereas the added term, Lj, 
is taken as the counter term that will be treated perturbatively. The mode expansion of the 
transverse vector potentials with the circularly polarized states are given by 

{ h+)kV(°+).(r)eS + &(-)kV^(°-).(r)4^] e^'-'' + hx.} , (40) 

(27r)2 ^2k{To) 

where the circular polarization vectors can be constructed out of the linear polarization 
vectors e^l and e^j as follows: 

4V = ^ + 4'.') . 4l' = {ie^l - 4'.') . (41) 

The corresponding mode equations obtained with respect to the unperturbative Lagrangian 
density then become 



CkV^%{r) = 0. (43) 

Thus, as compared with the mode equations of free photons ( l33l) and ( IMI) . the involved 
interaction gives the additional time dependent frequency terms in the mode equations that 
depend on the dynamics of the expectation value of the neutral pion field. In particular, 
when the mean field undergoes the oscillations around the equilibrium value to be determined 
dynamically, it is known that the solutions will exhibit the features of the unstable bands 
and the growth of the fiuctuating modes. The growth of the modes in the unstable bands 
translates into the profuse particle production. Thus, the photon production mechanism is 
that of parametric amplification. Here the background expansion will bring to us the extra 
damping effect. This is an extension of our previous studies on parametric amplification 
of the photon production to the case with the longitudinal background expansion. A novel 
phenomenon can be observed. We also observe the polarization asymmetry in the produced 
circularly polarized photons as a result of the pseudo-scalar nature of the coupling that is of 
interest to study. It is conceivable that after the heavy ions collision, the produced photons 
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are mainly having large momentum along the beam direction and then their kj^/k can be 
perturbatively small. The perturbation expansion will be adopted for this small parameter 
k']_/k'^. The interaction Hamiltonain density Hj can be constructed from the associated 
Lagrangian density L/ given by 



i^/k(r) 



1 ^ fc^ 



&(+)kV(°+),(r)-6(_)kl^(%(r) 



X fk ^ -kl . 



(44) 



Besides the interaction Hamilton depends on the amount of polarization asymmetry of the 
produced photons obtained from the mode equations P2l) and (113!) that certainly give the 
next order results. 

/^From the exact expressions of the electric and magnetic fields in Eqs. fl26l) and fl27|) . to 
order of k^/k"^ they can be obtained perturbatively as 



Kir) 



El^iT) + 5EliT) + 2t dr'eiT-r' 



BI^{t) + 2i I dT'e{T-T') BI^{t), H,^{r') +0{k\/k'' 



0{ki/k'), (45) 
(46) 



where 



— V-9{ 

2kiToy ^ 



-)]},(47) 



2k{ro) 



2k{To) 



V^k { [ (6(+)kV(°+).(r) + &|+)„k^(°4.(r)) eW^] - [ ( 



— 



(48) 
)] } (49) 



Then, the backreaction effects from the photon fields to the evolution of the mean field of 
the neutral pion in Eq. f lTT]) can be expressed by the electric and magnetic fields that in turn 
can be obtained perturbatively with all above results. They are 



EUr),B^Ur)\) + {{SE^o^.ir),B^Ur 



+2i / dr'eir-T' 



El 



fr), iJr 



/-k(T" 



+2i J dr' eir -t'){{ Bl^{r) , [Eo-k^(r), Hi.^{t') ] } ) + 0{ki/k') , 
9.8 {F,,F^p){t) = ( { B^^ir) , B^U-r) } ) 
Bi^^ir) i?o-k^(r) )+2ijdr'9{r- r') ( { Bl^{r) , [i?o-k^(r), i//„k(r ) ] } + 0(A;l/fc^ 



d_ 
d^ 



BI^{t) , B^Ur 
10 



(50) 



where the above correlation functions can be obtained from Eqs. (H7|) . ( HHj) . and (jUj), ex- 
pressed in terms of the mode functions, as follows: 



({^ok(r),5-k.(r)}) 



1 , d 
—n — r Q k — 

2k{To) dr 



[i?o-k.(r),if/_k(r')]}) 



2A;2(ro)^ Akr^ 

^{v^-)u{r)V^-W)-c.c.) - [(-)^(+)]} 

^{v^-)u{r)V^M^')-c.c.) - [(-)^(+)]} 



k' 



X 



[^(°).(r)V^&(r') 



c. c. 



(51) 



The curly bracket in above means the anticommutator. The results from the terms of order 
k^/k"^ give an estimate on the effects from their higher order contributions. 

According to the Bjorken's scenario the boost invariant hydrodynamical flow might be 
created after the heavy ion collisions. The corresponding initial vacuum state may already 
contain particles with respect to their asymptotical states to be observed in the detector j^. 
Besides the relaxation of the DCC can give an additional effect to produce the photons. Let 
us now consider the initial state at time tq given by the adiabatic modes in the comoving 
spacetime time that will be specified below. The corresponding initial particle number 
density can be expressed as 



(A4)(ro) 



2A;(ro) 



1 . 



(52) 



The expectation value of the number operator with respect to an initial vacuum state evolves 
in time and has the following form: 

( A4 ) (r) = [ ( Elir) E^Ur) ) + ( Bl{r) B^^r) ) ] - 1 , (53) 

that is, to order kj_/k'^, 

(A4)(r) = (Ei^^ir) Eo^Ur)) + {BUr) Bo^Ur)) + {{EUr) , SE^^.ir)}) 
+2i J dr'eir - t') {{El,^{r) ,[E,^Ur). ^7-k(r')] } ) 

+2z j dr'e{r -r'){{ Bl^{r) , [fio-k.(r), i//-k(r')] } ) + 0(^1/^^154) 
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Except for the correlation functions that we have found above, the correlation functions 
between the electric fields are obtained as 



({^ok(r), [^ok.(r),if,.(r')]}) = ^^^t^ E [ + c. c. 



X 



[V(°).(r)V^(°)*,(r')-c.c.] , (55) 



Thus the phase space number density in the comoving frame is given by [10|, lU 

dr](Px±dkrid'^k± ^■^'^'^ ( ) 

We now need to relate this quantity to the invariant spectra of the produced photons mea- 
sured in the laboratory frame. From the coordinate transformations (E]), we find that 

kr = p± COsh(77 — w) , 

k^ = —p±T sinh(?7 — w) , 

kx = Px, (57) 

where /c^ = {k, k) is the photon four-momentum in the comoving frame and = {p, p) is 
that measured in the laboratory frame. We here also introduce the outgoing photon particle 
rapidity defined by the four momentum in the center of the mass coordinate system given 
by 

P^J. = {p± cosh w , p± , p± sinh w) . (58) 

Hence, from the momentum transformation laws, one can change the momentum variables 
from kfj, into in the spectral number density given by 



dN dN 

P 



I dw dz drid'^x i dk„d'^k i ^ ^ 



dN 



d^p dwd'^k± J dw dz d-qd'^x^dkr^d'^k 

Thus, carrying out the Jacobian, the end result of the invariant spectra in the laboratory 
frame is obtained as 



— - -A f dk ^iV .gQ^ 

^ d^p dwd'^k± ^ J ^ d7]d'^x±dkrjd'^k± ' 

where A± is the transverse dimension of the effective transverse size of the colliding ions. 
This quantity is independent of as a consequence of the assumed boost invariance. 

Here we will postpone the full dynamical study on the photon production to our future 
work and consider free photons instead. With respect to an initial vacuum state of the 
photons given by the adiabatic modes, the particle number can also be measured by their 
asymptotical states to be observed in the detector under an expanding background. Exact 
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analytical solutions to the mode equations Eqs.( l37|) in the free photon case can allow to give 
an order-of-magnitude estimate on this effect, and they also provide an interesting example 
to illustrate the above formalism. To do so, we express the expectation value of the number 
operator as follows: 



T r • 



+ A;2(%),.(r)5(V),(r) + 5(.),(r)5(*_),(r))] -1. (61) 

The solutions S'(±) k) to the equations (!37l) can be found analytically in below 17| : 

Si±)kir) = q±)ifce-^i7«(fcxr) + C(±)2fce-^i7£)(fcxr) . (62) 

Here H^l^ and H^'^^ denote the Hankel functions of the first and second kind with purely 
imaginary order, u = /c^. Notice that the two solutions are complex conjugates of each other. 
Thus, the above constants Ci and C2 are to be determined by the initial conditions. They 
can be obtained from Eqs. (1351) and fl36l) in which the initial value of the mode functions 
V(i)fc and V(2)fc are specified by the adiabatic modes as 

-ifc{To)ro -ik{To)To I U2 \ 

V,„(.„) = = (-*(ro) -j^J , (63) 



g-ifc(To)ro _ g-ifc(ro)To 

^(2)fc(ro) = ^— , V^2)k{ro) = ^—{-ik{To)), (64) 



consistent with the Wranskian conditions (12 9 p and (130 p . Fig. 1 displays the result of the 
spectra of the photons per unit effective transverse area of the collisions associated with 
its initial state measured with the reference of the asymptotical states. Typically, less 
one particle per unit effective transverse area is found [l6|. It would be of interest in both 
theoretically and experimentally to learn how the photon production can be further amplified 
from the relaxation of a DCC under a longitudinal background expansion. 

These derived dynamical equations display several important physical processes: i) non- 
linear relaxation of the mean field within a DCC, ii) particle creation due to parametric 
amplification of the fluctuations driven by the time dependent mean field, iii) the enhance- 
ment from the boost effect under a longitudinal hydrodynamical expansion. This must be 
studied self-consistently from the coupled set of Eqs. f ITT]) . ([HD, and fISUl) by inserting the 
mode functions in Eqs. fITH]) . ( W2\i . and (143 p . The invariant spectra of the observed photon 
production can be obtained from Eqs. (1541) . (!55l) . and (!60|) . In particular, the expressions 
of corrections give an estimate on the ignored terms with higher order k^/k"^. Although 
we do not intend to offer a numerical analysis of the resulting equations in this paper, the 
nonequilibrium equation of motion has the potential for providing an enhancement of the 
photon production. This enhancement could lead to an experimentally observable signal in 
the direct photon measurements of heavy-ion collisions that can be a potential test of the 
formation of disoriented chiral condensates. We postpone to a forthcoming paper the full 
numerical study of these equations and an assessment of the potential phenomenological 
impact of the nonequilibrium dynamics. 

In conclusion, what we have done in this work is to propose the quantization scheme 
on electromagnetic potentials under a longitudinal background expansion, which suits the 
nonequilibrium initial-valued problem. This is a necessary first step toward the study of 
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1 2 3 4 5 



FIG. 1: The respective spectra of the photons per unit effective transverse area of the cohisions 
in its initial state measured by the asymptotical states as a function of the transverse momentum 
( in units of tq) at times ( in units of tq ) r = 2 ( long-dashed), 5 ( short-dashed), 10 (solid). 



photon production through the non-equihbrium relaxation of a disoriented chiral condensate 
formed in ultra-relativistic heavy-ion colhsions. The next calculation to do would be a 
numerical study of the above equations with an eye on obtaining the spectra of the produced 
photons that can be compared with thermal photons from quark-gluon plasma and hadronic 



matter as done in Ref. This would have fascinating phenomenological consequences. 
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